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Abstract 

We present a model based on the flavour group S3 x Z3 x Zq to explain the main features 

of fermion masses and mixing. In particular, in the neutrino sector the breaking of 

O y the S3 symmetry is responsible for a naturally small r = Am^/Am^ and suitable 



next-to-leading order corrections bring #13 at the level of ~ 0.13, fully compatible with 



D ■ the recent Daya Bay result. In the quark sector, the model accommodates the different 



mass hierarchies in the up and down quark sectors as well as the Cabibbo angle and 
I Vcb I ( or I V u b I , depending on the charge assignment of the right-handed b-quark) in the 
correct range. 
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1 Introduction 

One of major challenges in the model building is to provide a framework where the 
hierarchies in the fermion masses and the values of the mixing angles can be explained in a 
natural way. Recent results in the neutrino sector [TJ have inspired many authors to discuss 
the possible origin of a relatively large reactor angle, as obtained first by the T2K and 
MINOS collaborations j2] and, more recently, by Double Chooz [3] and Daya Bay jl]. In 
particular, the latter has given a 5.2a evidence of #13 ^ 0: 

sin 2 26 13 = 0.092 ±0.016 ±0.005, (1) 

probably solving the longstanding question of the magnitude of #13 (see [5], [6], [7], [8] 
and [9] for a discussion of the impact of such a result for lepton flavour mixing and leptonic 
CP violation). Several models based on the breaking of additional symmetries beyond the 
Standard Model have been proposed to explain this and other peculiarities emerged from 
neutrino oscillation experiments, as the large solar and atmospheric mixing angles as well 
as a small solar-to-atmospheric mass differences ratio r. In an unified view of fermion mass 
hierarchies and flavour mixing, the quark sector must also be taken into account and this 
poses the question of how to reconcile the numerous differences in the spectra and mixing 
of quarks compared to leptons, the most relevant ones being a different structure of the 
CKM mixing matrix and the strong mass hierarchy in the up-quark sector. Although many 
attempts have been done in this direction, the question of explaining the values of fermion 
masses and mixing (one of the aspects of the flavour problem) is still an open issue and 
deserves further studies. Discrete symmetries [TU] have been invoked as a powerful tool to 
solve/mitigate some of these problems and, among several choices, the non-abelian discrete 
group S3 [11] has also been investigated. Remarkably, models based on S3 describing both 
quarks and leptons (with very few assumptions on the scalar and Yukawa sectors of the 
theory) are not very common. Interesting attempts in this direction have been done in the 
context of GUT theories [TJ] and in a non-unified approach [TJ]- [14] . In this paper we want 
to contribute to this discussion, presenting a SUSY see-saw model based on the flavour group 
S3 x Z 3 x Zq. Identifying e as the small S3 breaking parameter (whose magnitude will be 
discussed later), the main features of our construction are: 

• the neutrino mass hierarchy can only be of normal type; 

• the solar-to-atmospheric mass differences ratio r is naturally small, at the level of 
0(e 2 ), without invoking any fine-tuning among the Yukawa parameters (Sectf2J); 

• after the inclusion of suitable next-to- leading order corrections (NLO), the vanishing 
leading order (LO) #13 receives a relatively large shift of 0(e) whereas #23 is only 
modified by terms of 0(e 2 ); the solar angle remains large, of 0(1) (SectsflJlHl); 

• in the quark sector, we are able to reproduce the correct order of magnitude of the 
mass hierarchies in both up and down sectors (Sect|H]); 

• the Cabibbo angle is predicted in the correct range as well as \V c b\ (or \V u b\, depending 
on the charge assignment of the right-handed b-quark); 

• the flavon alignments needed to get the previous results are obtained from the mini- 
mization of the superpotential allowed by the symmetry of the model (SectfJJ). 



A bit more problematic is the explanation of the muon-to-tau mass ratio and the value of 
| V u b | (or |Vd,|), in any case naturally smaller than the Cabibbo angle. Other interesting 
features of our model are the prediction for the effective mass \m ee \, quite small as usual for 
neutrinos with normal mass ordering, and the fact that the Tri-bimaximal mixing structure 
(TBM) [15] can be easily obtained using only one additional relation among the Yukawas 
(not dictated by the symmetry of the theory). 

The model is formulated in the framework of the SUSY see-saw mechanism. At LO, the 
contribution to # 12 only comes from the neutrino sector whereas #23 is completely generated 
by the charged lepton sector. At this order, the reactor angle is vanishing. To produce such 
a pattern, the mass matrices mim\ and m u must be block-diagonal in the (23) and (12) 
sectors, respectively. In particular, having a block-diagonal light neutrino mass matrix helps 
in giving different order of magnitude to the mass eigenstates and then in obtaining a small 
r . 

We use a non- conventional assignment of the left-handed doublets L of the second and 
third generations in a 2 representation of S3 as follows: 

;)< < 2 » 

whereas for the right-handed doublet we assume: 

£ ) • (3) 

Electron fields are assigned to the singlet 1. In the neutrino sector, we introduce three 
right-handed fields: the first two generations are grouped in the 2 representation: 

% ) • 

whereas z/| belongs to a I2. We use the following representation of the S3 group: 

S 2 = T 3 = (ST) 2 = 1 . (5) 

In the so-called "T-diagonal" basis we have: 

*-(;;) '-(o£). « 6 » 

with <jj = e* 2 / 37r . The tensor products involving pseudo-singlets are given by I2 x 1 2 = 1 and 
I2 x 2 = 2 while the product of two doublets is 2x2 = 2 + 1 + 1 2 which, in terms of the 
components of the two doublets A = (01, «2) T and B = (b±, b 2 ) T , are as follows: 

ai b 2 + a 2 b x e 1 a x b 2 - a 2 b x e 1 2 ( ^ 2 J G 2 . (7) 

To ensure the breaking of S3 along appropriate directions in the flavour space, we need 
two doublet flavons and £ with the following vevs (derived from the minimization of the 
superpotential, see Sectj3]): 

(*> = «*(!) (0 = «c(J); (8) 

two other scalar singlets x ~ 1 and x' ~ I2 are also needed to guarantee appropriate non- 
vanishing entries in the fermion mass matrices; both singlets have non vanishing vevs: 

(X) = v x (X') = v' x . (9) 

The particle content and the transformation properties of leptons, electroweak Higgs doublets 
and flavons under S3 x Z3 x Z s are summarized in TabCEJ 
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Table 1: Transformation properties of leptons, electroweak Higgs doublets and flavons under 
S 3 x Z 3 x Z 6 and U(1) R . 

2 Leptons: leading order 

2.1 Charged leptons 

The leading order lagrangian in the left-right basis reads: 



(p v 

C d = x 1 LL c h d —+x 2 LL c h d — . 



(10) 



At this stage, there are no contributions to the electron entries (since x' has been assigned 
to the I2 representation) and the electron mass is vanishing. In the /i — r subsector the mass 
matrix is as follows: 

/ X 2 V X XiV<j, 
XlV<j, x 2 v x 

where v d = (hd)- The phases of the Yukawas x\ )2 can be rotated away, with no loss of 
generality. The charged lepton masses and the mixing matrix are then: 
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(12) 



(13) 



it has been 
, however, 



The question is how to realize a naturally small muon-to-tau mass ratio. In 
suggested to introduce a further symmetry under which x± = —x 2 and v^ == v 
as pointed out in [14] , this additional symmetry does not commute with S3 and the whole 
symmetry group is larger than the latter. Here we observe that for real flavon vevs (as 
assumed here) any couple of real (xi,x 2 ) determines a larger and a smaller eigenvalue, so a 
mass splitting is quite natural in this model. If we want to correctly reproduce the mu-to-tau 
mass ratio we have to fine-tune the Yukawa parameters; in particular, assuming m^/m T ~ A 2 
(A being the Cabibbo angle) we get: 



•'"1 



-x 2 (l + 2\ 2 ) + 0(\ A 



(14) 



(so that, as expected, the two Yukawas must be opposite in sign and almost equal in mag- 
nitude). We notice that, at this level, we do not need to specify a value of the ratios of the 
flavon vevs over the cut-off A. 



2.2 Neutrinos 

The first useful terms to generate a Dirac mass matrix are: 

x' 

C D = aLv c h u + bL e vl~h u — , 



(15) 



so that: 
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(16) 



where v,, 



and: 



(h u ). For the Majorana mass matrix we have: 
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The first contribution to the heavy masses arise from the previous matrix with v' 
are: 




(17) 



(18) 



and 



(19) 



The degeneracy is then lifted by taking into account the corrections from the v' terms. The 
light neutrino mass matrix is obtained from the see-saw formula: 



m,. 



-m D mJmT 



"M " b D i (20) 

to make the following analytical evaluations more readable, we fix a = 1 (but it will be 
considered as a free parameter in the numerical evaluations to follow) and assume an equal 
order of magnitude of the flavon vevs to the cut-off ratio M- = ^ = e, that we take as a 
small parameter. At first order in e we get: 
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(21) 



All phases can be reabsorbed into a redefinition of the right-handed neutrino fields, so we 
deal with real parameters. We see that, at this order, we get a block- diagonal form of the 
light mass matrix, with the contribution to the (12) sector of 0{e) and the (33) element 
larger. The matrix m u is compatible with normal hierarchy only and the light masses (still 
at the first order in e) are then given by: 
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so that the solar and atmospheric mass differences read: 
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(23) 
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Taking a = 1, Am^, = 2.5 x 10" 3 eV 2 and v u = 100 GeV, we estimate d$~2x 10 14 GeV, 
which is a common order of magnitude for the heavy neutrino masses. It is now easy to 
derive the ratio r: 

_ AmL _ c^(4& 2 /3 2 + ^/ 2 ^ 2 

~ Ami" 0* (24) 

which, for 0(1) parameters, is naturally suppressed by e 2 . The value of e cannot be precisely 
determined at this stage; to recover the experimental value r ~ 1/30, e should be as small 
as the Cabibbo angle but a small numerical enhancement (suppression) of the coefficient 
- — - — T^r 1 — - — can bring it to smaller (larger) values. The matrix in eg. ( 12 1 \ determines a 
non-vanishing solar angle: a straightforward computation gives: 

^ 3 ^ = 1 - WP + 8W+8 ' (25) 

so that it is generically of 0(1). Summarizing, the whole neutrino mixing matrix at LO 
Upmns — UJ U u has a vanishing reactor angle, maximal # 23 and large 0(1) solar angle. 

It is interesting to observe that, to reproduce the TBM value tan 2 # 12 = 1/2, one simply 
needs 



5 = ±b/3/V2, (26) 

which is still a number of 0(1); in this respect, this external condition does not appear to 
be completely unnatural since it does not require any strong hierarchy among the model 
parameters. Another interesting possibility to get the TBM matrix from the S3 symmetry 
has been proposed in [17] where, however, after the inclusion of the charged lepton corrections 
one of the two allowed S3 invariant Yukawa couplings must be switched off by hand. Notice 
that a maximal value for 6* 12 in our model can only be obtained if 5 = 0, which implies a 
vanishing Am 2 



sol' 



3 Flavon alignment 

The structure of the flavon vevs can be obtained minimizing the scalar superpotential 
in the limit of exact SUSY [18]. Within this approach, a continuous U(1)r symmetry is 
introduced, under which matter fields have R = 1, while Higgses and flavon fields have 
R = 0. Such a symmetry will be eventually broken down to the R-parity by small SUSY 
breaking effects which can be neglected in the first approximation. Since the superpotential 
must have R = 2, we need to introduce two additional scalar fields, a doublet Ao and a 
singlet £o 5 with R = 2. Within these assumptions, the relevant part of the scalar potential 
of the model is given by the F-terms with: 



i 

In the following, we parametrize the vevs as: 
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difi 



(27) 



(0) = v, ( £ j (0 = v, ( I J ( X ) = v x (x') = v' x , (28) 

where 0i j2 and £i j2 are adimensional quantities. At the leading order we have: 

w d = gi ^e + 92^o(f> 2 + 93^o(l)X- (29) 



The condition for the minima are: 
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92 V f 02 + 93 v<t, (f) 1 v x = (30) 



66 = 0. 



6 
The set of equations admit the solution: 

0i = 02 v x = v$ 6 = 0. (31) 

93 

The relation among v^ and v^ allows us to assume a common order of magnitude for these 
vevs; on the other hand, the choice 6 = 0, 6 7^ is not restrictive since the other one, with 
6 7^ 0,6 = [13] , can be obtained acting with the generator S on it. At this stage, the 
flavon x' does not appear in the superpotential. However, the first corrections to Wd involve 
X' A and read: 

5w d = g[Z ^. (32) 

Since it appears with the driving field 6> this term can modify the vev of the £ flavon; 
assuming a perturbed structure like: 

(0 = vt ( j « g ) , (33) 

the minimizing equation (at the first order in the perturbations £[ 2 ) gives: 

/ 2 \ 2 

Zgivlti + til^f) =0. (34) 

Dividing this equation by 1/A 2 and assuming, as usual, that v^/A ~ v x r/A ~ s, we can 
estimate 6 ~ £<2 an d a l so: 

■fc'-ff- (35) 

The perturbation £ 2 remains unspecified and we put it to zero. Summarizing: 

(0 =vJ ^f ), v' x ~ « c , (36) 

with 6 a coefficient of 0(1). For the other two flavons and x> the first useful corrections 
arise at the level of five-flavon insertion. The corrective terms are then of relative 0(e 3 ) with 
respect to the leading order results and will be neglected in the following. Then, the NLO 
corrections to the mass matrices will be computed using the vev structures for and x as 
given in eqs.®-©. 



4 Next to leading order corrections 

It is important to check that the previous results on the mixing angles are not destroyed 
once the corrections to the lagrangians are taken into account. 



4.1 Charged leptons 

The most relevant corrections come from the term: 



5C d = x 3 L e L c h d 



A 2 



(37) 



which modifies the mass matrix as follows: 



x 3 e 2 —x 3 e 2 






(38) 



Still working with real parameters for simplicity, we see that the expression of the masses 
are not modified whereas the mixing matrix is now given by: 
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(39) 



Then, there are new contributions to the leptonic 0± 2 and 8 13 but not to the atmospheric 
angle, since the 2 — 3 sub-block still remains unchanged at this order. The electron mass is 
not generated at this order but only at the level of 5-flavon insertion, from the following two 
(non-vanishing) operators: 

L e L c e h d {e<P 2 xA 2 <Pxx) , 

which give a contribution of 0(e 5 ); reminding that m T ~ e, the electron-to-tau mass ratio 
is of 0(e A ). Considering that this ratio should be some units of 10~ 4 , we deduce that e 
cannot be too different from ps 0.13. It is important to stress, however, that this value of 
e is merely indicative, as it can be enhanced or suppressed by a proper arrangement of the 
Yuwaka couplings. We will study more in details the magnitude of the breaking parameter 
in Sectj5l As it will be clear in the next section, Uf will provide the largest corrections to 
the LO results. 



4.2 Neutrinos 

At the next level of approximation, the light neutrino mass matrix can be still evaluated 
using eq.( l20~]) but expanded now at 0(e 2 ) with flavon alignments as discussed in Sectj3j 
Here we want to comment that the filling of some of the vanishing entries in eq. (fT6|) and ( TT8|) 
requires multiple-flavon insertions; for example, the elements (22), (23), (31) and (33) in mo 
are generated by operators with three flavons (like Li> c <j) 3 h u ) and are then of 0(e 3 ). All 
next-to-leading order Majorana operators are also of 0(e 3 ) since they contain at least four 
flavons. Then, the main contributions to the LO vanishing matrix elements in eq. ffTB"]) come 



from the vev shift in eg.( !36p and are of 0(e 2 ). From these considerations is not difficult to 
understand that the NLO light neutrino mass matrix is given by (still using a = 1): 
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(40) 



The expression of the neutrino mixing matrix is quite cumbersome; we have found that all 
eigenvalues and eigenvectors are corrected by (intricate) 0{e 2 ) terms. As previously stated, 
the charged lepton rotation Ug gives the main corrections to the LO results for #i 2 and #13, 
whereas the atmospheric angle only receives 0(e 2 ) corrections from the neutrino sector. The 
final results for the mixing angles are then: 
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(41) 
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As we can see, it is not difficult to reconcile our results with the experimental data, barring 
accidental cancellations; in fact, the NLO have preserved many good features of the LO 
results (large solar and atmospheric mixings) while producing a relatively large shift for 
6 13 ~ 0(e). 



4.3 Effective mass terms and 
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The previous results are not substantially modified by effective Weinberg operators. Up 
to four-flavon insertion, the lagrangian is as follows: 
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where Aln is the lepton number breaking scale. The contributions to the neutrino mass 
matrix is then: 
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(43) 



so we still have a block-diagonal form. In practise, the relevant contribution to m u is given 
by the first operator in C w , which fills the (11) vanishing entry in eg. (HOI) with a term of 
0(v 2 /Aln)- In the case Aln ~ A, this term also contributes to eg. (1221) and eg.(|23l). changing 
the coefficients in front of the e parameter but not their order of magnitude. In the case 
Aln S> A, this term is negligible and the effective operators do not play any relevant role in 
the determination of the neutrino masses and mixings. 

To evaluate the prediction of our model for the effective mass m ee , we work in the basis 
where the charged leptons are diagonal and extract the (11) matrix element of the rotated 
neutrino mass matrix, that is: 
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We get: 



\m P 



sin 2 9 



13 



[2ab(xi - x 2 ) + /3x 3 } 



a (3 x 3 



(45) 



is small, at the level of 9f 3 , although 



As expected in models for the normal hierarchy, \m ee 

a clear correlation with the reactor angle is lacking because of the 0(1) coefficient. 



5 Numerical analysis of the lepton sector 

The main purpose of this section is to analyze in detail the implication of our model 
for the lepton masses and mixing. In doing that, we use the NLO charged lepton and 
neutrino mass matrices as given in eqs. (TTT]) - (l38|) and eg. (140]) and perform a MonteCarlo 
simulation extracting complex lagrangian parameters with absolute values in the interval 
[1/2,2] whereas the small S 3 breaking parameter is taken randomly in [A 2 , A]. To study 
more in detail the magnitude of e, we imposed relaxed bounds on the charged lepton mass 
ratiocl: 



Tfl 

2 x 1(T 3 < — < 9 x 10" 3 



1-2 



Tfl 

2 x 10~ z < — £ < 9 x 10" 



(46) 



we also impose the constraints on the neutrino mass differences: 
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; eV 2 , (47) 



derived from the second reference in [I], so that the value of r is automatically reproduced. 
The main results of such an analysis are presented in FigsJTJEJ In the former, we are interested 
to the dependence of sin 2 #12 and sin 2 #23 on the breaking scale e. In these plots, the horizontal 
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Figure 1: Mixing angles sin 2 6* 12 (left panel) and sin 2 6*23 (right panel) as a function of the 
breaking scale e. The horizontal gray bands are the regions excluded by the experimental data 
at 3a, as obtained from the second reference in FT]/. 



gray bands are the regions excluded by the experimental data at 3a, as obtained from the 
second reference in [T]. We clearly see that e ~ 0(\ 2 ) are almost excluded because too small 
to fill the relations in eg. (146 p . For every mixing angle, the bulk of the selected points is 
around e ~ 0.14, not really different from the estimate we gave in Sect J4. II The analytical 



2 More restricted bounds would only select a narrower range of e's. 



results anticipated in eq. (l4"Tj) are also confirmed; in particular, we see that the solar angle 
(left plot) is mainly undetermined but a large fraction of the points fall into the allowed 3a 
range, as a consequence of the fact that the TBM approximation is, in some sense, contained 
in the model via the simple relation in eg. (126]) . For the atmospheric angle (right plot), the 
majority of the points are well inside the 3a range, showing also a tendency to a largest 
spread for e > 0.14, as indicated in eq.( )4TT) . Similar considerations could also be drawn 
for the dependence of the reactor angle on e; however, insted of presenting a scatter plot, 
we prefer to compare the sin 2 2# 13 distribution as obtained from our numerical simulation 
directly with the Daya Bay result. This is shown in Figj2j where the 3a bounds derived 
from eq.([T]) are enclosed in the solid vertical lines, whereas the dashed line is the best fit 
point. We can appreciate that, although the distribution is quite broad, the largest density 
of sin 2 2# 13 extractions is just inside the allowed region. We stress that this result has been 
obtained with an ab-initio simulation of the charged lepton and neutrino mass matrices with 
no other constraint than those given in eqs. (j46ti47|) . Notice that, given the large number of 
0(1) parameters, no definite predictions for the CP phase can be drawn. 



sin 2 (0, 3 ) 



Figure 2: Distribution of the valus of sin 2 29 13 as obtained from our model; the vertical solid 
lines enclose the 3a range derived from eq. (QP whereas the dashed line is the corresponding 
best fit point. 



6 The quark sector 

The S3 x Z3 x Zq symmetry provides a good description of the quark sector also; we 
use the same flavon fields and alignments described in the previous sections. The first two 
families of left-handed quarks are assigned to a 2 representation of S3 whereas all other fields 
belong to singles 1 or 1 2 in the case of b c (see Tabf2]). The lagrangian in the up-quark sector, 
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Table 2: Transformation properties of quarks under S3 x Z3 x Zq. 
including all relevant operators to generate non-vanishing entries in the mass matrix, reads 
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as follows: 
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where we have indicated the S3 contractions when necessary. The related mass matrix is as 
follows: 
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2 c u £' e 4 b u e 3 a u 



First of all, we observe that the mass hierarchy is well reproduced for the same value of e as 
deduced from the lepton sector; in particular, we have: 

m u = v u (f u + f' u ) e 5 m c = v u e u e 3 m t = v u a u . (50) 

Then, taking real Yukawas for simplicity, the matrix U up diagonalizing m up m u is given by: 



1 


e' 


du £ 3 


du £ 3 


1 
du £ 3 


_du £ 3 
1 



U up = -i e 1 -^e 3 • (51) 



This matrix goes into the desired direction: at LO, U up is the identity and at the NLO there 
exist a well defined hierarchy among the (12) element, not far away from the value of the 
Cabibbo angle, and the other off-diagonal matrix elements, which contribute to V c t and V u b. 
From the charge assignment in TabfSJ we see that the operators in the down sector involving 
d c and s c are the same as those with u c and c c , respectively (with the obvious replacement 
h u — > hd), whereas all operators including b c do not have a corresponding in C up ; then we 
have: 

r h ucX' . , , c i 2( i> . , id 2 . , n , , c (0 2 )i^ 

£down = a d q 3L h d b — + b d q 3L h d s —j- + c d q 3L h d d —^ + d d Q L h d b 

+d' d Ql h d b c ^^ + d' d Q L h d b c ^ + d" d ' Q L h d b c ^ (52) 

+e d Q L h d s — — + e d Q L h d s -— + f d Q L h d d — 

with mass matrix as: 

(f<i + f'd)/ ~e' d e 3 -(2d d + d' d )e 2 
mdown = v d e ~(2f d + f' d ) e A e d e 2 (d' d + d d ) e 2 (53) 

2 c d g e 3 b d e 2 a d 
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Again, the mass ratios are well reproduced, since: 



m d = v d (f d + f d )e 5 m s = v d e d e 3 m b = v d a d e 




(54) 


Also, the bottom-to-top mass ratio is given by 




— = — £ (tan/3) , 
m t \a u J 


(55) 


implying tan/3 ~ 5 for Yukawas of 0(1). The matrix U down diagonalizing m down 


m Lu,n 1S 


given by: 




/ 1 t £ (2^ + ^)e 2 ) 






U *=\ -%z X -K + ^> 2 




(56) 


\-(2d d + d' d )e 2 «, + <£V 1 J 







Two comments are in order: on the one-hand, the (12) element is still of the correct order 
of magnitude to explain, in combination with the result from the up sector, the value of the 
Cabibbo angle. On the other hand, all other off diagonal entries are smaller that (12) but 
the (13) element is a bit larger than the required values to fit V ub . In fact, these elements in 
eq. (!56j) are the dominant contributions to V cb and V ub in the CKM, since the corresponding 
matrix elements in eq. (j3T|) are generally smaller. In fact, we get: 



(57) 



(---) e V ub = (2d d + d d )e 2 



Vcd — —Vus Vtd — —Vub Vts — —Vcb ■ 

Obviously, the last equalities are lifted once the 0(e 3 ) terms are taken into account. To 
corroborate the previous considerations, we perform a numerical simulation of the up and 
down mass matrices of eqs. (149p and ( 1531) extracting, as we did for the lepton sector, complex 
Yukawa parameters with absolute values in the interval [1/2, 2]. We fixed the breaking scale 
e = 0.15, as suggested by the same procedure done in the lepton sector. We also impose the 
following constraints on the relevant mass ratios [T9] : 



(58) 



Tfl Tfl 

5.40 x 10~ 6 < — < 1.18 x 10~ 5 3.60 x 10~ 3 < — < 4.10 x 10" 
m t m t 

Tfl si TH 

7.90 x 10" 4 < — < 1.34 x 10" 3 1.60 x 10" 2 < — < 2.90 x 10" 



m b m b 

Our results for every CKM matrix elements are shown in figJSl where we plot the absolute 
values of the nine distributions of the Vij entries. The diagonal entries are displayed in linear 
scale whereas we adopt a log scale for the off-diagonal elements; for them, we also showed 
the best fit values [20] with a solid vertical line. As we can see, it is very easy in our model 
to reproduce, with no fine-tuning, the correct values of V us , V cd as well as V cb , V ts while a 
small discrepancy remains with the best fit of V ub and V td , as anticipated above. To make 
also these matrix elements fully compatible with the data, we need a moderate cancellation 
among the Yukawas d' d and d d , see eq. (156]) . We do not present any plot related to the CP 
phase since the large number of (9(1) parameters does not allow any definite prediction. As a 
last comment, we observe that taking a charge assignment for the field b c as the one adopted 
for t c , we would get a structure of U d similar to U up as given is eq. floTl) : this automatically 
would imply V ub ~ V cb ~ V td ~ V ts ~ e 3 , so an enhancement should be invoked to reproduce 
the values of V ub and V td . In this case, the bottom-to-tau mass ratio would be completely 
explained by large tan/3. 
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Figure 3: Distributions of the CKM matrix elements as obtained in our model. The best fit 
values WD^ of the off-diagonal matrix elements are shown with a solid vertical line. 



7 Conclusions 

The Daya Bay Collaboration has released the measurement of the reactor angle #13, 
showing a 5.2a discrepancy from zero. From the model building point of view, neutrino 
mass textures predicting a vanishing #13 at leading order seem to be less appealing, unless 
large corrections bring the reactor angle to values compatible with the recent results, without 
destroying the predictions for the other mixing parameters. In this paper we have presented 
a see-saw SUSY model for fermion masses and mixing based on the non-abelian discrete 
symmetry 53, whose main result is the prediction of a large #13 ~ 0.13, fully compatible 
with the Daya Bay claim of eq.fCQ). Other remarkable features of our construction are: 

• in the lepton sector, the spectrum is of normal type, with Q^i an d 6*23 compatible with 
their experimental allowed ranges; 

• in the quark sector, we obtained a good description of the relevant mass ratios and 
the absolute values of all the CKM matrix elements (including the Cabibbo angle) but 
V u b, Vtd, for which we need a moderate fine-tuning among the Yukawas defining these 
matrix elements; 

• the flavon alignments needed to get the previous results are natural minima of the 
superpotential in the SUSY limit. 
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